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Abstract. The irreversible adsorption of polymers to a two-dimensional 
solid surface is studied. An operator formalism is introduced for chemi- 
sorption from a polydisperse solution of polymers which transforms the 
analysis of the adsorption process to a set of combinatorial problems on 
a two-dimensional lattice. The time evolution of the number of polymers 
attached and the surface area covered are calculated via a series expansion. 
The dependence of the final coverage on the parameters of the model (i.e. 
the parameters of the distribution of polymer lengths in the solution) is 
studied. Various methods for accelerating the convergence of the resulting 
infinite series are considered. To demonstrate the accuracy of the truncated 
series approach, the series expansion results are compared with the results 
of stochastic simulation. 



1 Introduction 

The adsorption of polymers to solid surfaces has wide technological and medical ap- 
plications |14l 1^. In this paper, we study chemisorption, i.e. the situation where co- 
valent surface-polymer bonds develop and adsorption is effectively irreversible on the 
experimental time scale |13j . Chemisorbing polymers have one or more reactive (bind- 
ing) groups along the polymer chain which can react with binding sites on the surface. 
Polymers with one reactive group at the end of the chain are called semitelechelic. A 
schematic diagram of the adsorption of a semitelechelic polymer is shown in Figure 
^a) where the binding sites are arranged into a rectangular mesh on the surface. An 
important parameter of the chemisorption process is the density of binding sites, or 
equivalently, the average distance between neigbouring sites, which is denoted by h in 
Figure^a). Denoting the hydrodynamic radius of the polymer by i?, we can distinguish 
three different scenarios. If <C i?, then the polymer layer created by chemisorption of 
the semitelechelic polymer will be a polymer brush after sufficiently long time jllll??!!^. 
In this case, one can simply assume that a polymer can attach anywhere on the sur- 
face for modelling purposes. In particular, one can use continuum random sequential 
adsorption to model the process The other extreme case is h ^ R where the final 
layer contains one attached polymer at each binding site. No steric shielding needs to 
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Figure 1: (a) j4 schematic diagram of the chemisorption of the semitelechelic polymer. 
(b) Situation from (a) transformed into two-dimensional lattice setting. 

be considered when modelling the process and the dynamics of adsorption is trivial from 
the mathematical point of view. The last important case is when h R. This is the 
regime studied in this paper. 

Chemisorption is often modelled as a random sequential adsorption (RSA) 01 . 
In a previous paper |3] we studied one-dimensional models of random sequential (irre- 
versible) adsorption. Our motivation was to understand the essential processes involved 
in pharmacological applications such as the polymer coating of viruses 'S'. The classical 
RSA model 4^ was generalized to study the effects of polydispersity of polymers in solu- 
tion, of partial overlapping of the adsorbed polymers, and the influence of reactions with 
the solvent on the adsorption process. Working in one dimension, we derived an integro- 
differential evolution equation for the adsorption process and we studied the asymptotic 
behaviour of the quantities of interest, namely the surface area covered and the number 
of molecules attached to the surface. We also presented applications of equation-free 
dynamics renormalization tools jlOj to study the asymptotically self-similar behaviour 
of the adsorption process. In [3] we used a continuum RSA model. The underlying 
assumption was that the polymer can effectively bind anywhere on the surface, i.e. we 
worked in the regime h <ti R. In reality, the reactive groups on the polymer can react 
only with the corresponding binding sites on the surface, which are primary amino- 
groups in the virus coating problem. Rough estimates from molecular models suggest 
that the average distance between primary amino-groups in the virus capsid is about a 
nanometre 0. However, it is difhcult to guess which of the amino-groups in the capsid 
are available for the reaction with the polymer, i.e. are accessible for polymers from 
solution. In particular, both the regimes h R and h ^ R can be justified in the 
virus coating problem. Other chemisorbing systems |14l ^ can be also used to motivate 
investigation of the borderline case /i ~ i?. 

Assuming h ^ R, we have to take the discrete nature of the binding sites into 
account. This means that lattice RSA modelling is more appropriate than continuum 
RSA modelling. In this paper we assume for simplicity that the binding sites lie on a 
rectangular mesh (see Figure ^1, with mesh points a distance h apart. Wc choose h — 1 
without loss of generality in what follows. Any polymer covers the binding site to which 
it is attached. Moreover, longer polymers also effectively cover neighbouring binding 
sites, as illustrated in Figure ^a). More precisely, an attached semitelechelic polymer 
covers a circle of a certain radius r which is centered at the binding site (meshpoint 
{i,j)). If r < 1, then the polymer effectively covers only the corresponding binding site 
{i,j). If 1 < r < \/2, then the polymer covers a small "cross" Xij where we define 



= {(l, j), {l + {l - {l,J - 1), {l,J + 



(I.l) 
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We call set of mesh points Xij the cross (or cross-polymer) centered at (z, j). If < 
r < 2, then the polymer covers a small "square" Sij defined by 

= {(«, j), (i + 1, j), (« - 1, j), («, j - 1), («, j + 1), 

{t + + 1), (z + 1, J - 1), (z - 1, J + 1), (z - l,j - 1)}. (1.2) 

We call set of mesh points Sij the square (or square-polymer) centered at (z,j). If 
2 < r, then the polymer covers at least 13 binding sites. To simplify the combinatorial 
complexity of the problem, we restrict our consideration to the case r < 2. In this case, 
we can formulate the chemisorption of polymers in terms of adsorption of points, crosses 
and squares to the two-dimensional lattice (see Figure^b)). We denote by a the fraction 
of polymers in the solution for which 1 < r < v^, so that a is the probability that a 
randomly chosen polymer in solution will adsorb as a cross. Similarly, we denote f3 the 
fraction of polymers in the solution for which < r < 2 so that /3 is the probability 
that a randomly chosen polymer in solution will adsorb as a square. In particular, 
we must have < a + f3 < 1 where 1 — a — /3 is the probability that a randomly 
chosen polymer in solution will adsorb as a point. We work with an M x AI mesh 
with periodic boundary conditions. Then our two-dimensional polydisperse random 
sequential adsorption (pRSA) algorithm can be stated as follows. 

pRSA algorithm: We consider the adsorption of points {{i,j)}, crosses Xij 
and squares Si_j to the two-dimensional rectangular M x M mesh. At each 
time step, we choose randomly a point {i,j) in the mesh. If the selected mesh 
point {i,j) is covered (occupied) by a point/cross/square already placed, the 
adsorption is rejected. If the mesh point {i,j) is vacant, then it is marked as 
occupied. Moreover, with probability a (resp. j3), all mesh points in the set 
Xij (resp. Sij) are marked as occupied. 

To simulate pRSA algorithm, we have to generate three random numbers at each time 
step. The first two of them are used for random selection of the lattice point where 
the reactive group of the adsorbed polymer is attempted to bind. The third random 
number r„, uniformly distributed in interval [0, 1], is used to determine the length of the 
adsorbed polymer. If r„ G [0, a), then the cross polymer is placed. If r„ £ [a,Q -|- /3), 
then the square-polymer is chosen. If r„ g 1], then the point-polymer is adsorbed. 

An illustrative numerical simulation of pRSA algorithm for a = 0.8 and /3 = 0.1 is shown 
in Figure 13 We start with an empty rectangular 100 x 100 mesh, i.e. M = 100. The 
mesh points covered by polymers are plotted at different times. 

Let us note that pRSA algorithm requires that the position {i,j) of the center of 
the adsorbed cross Xij (resp. square iS^^) is vacant. On the other hand, the "tails" of 
crosses/squares can overlap. Here, the center of the cross (resp. square) describes the 
reactive group which is covalently bound to the surface. The remaining four (resp. eight) 
points of the cross (resp. square) describe the polymer tails which sterically shield the 
neighbourhood of the adsorbed polymer. In our algorithm, binding of a larger polymer 
prevents binding (of the center) of another polymer in the neighbourhood of the center 
of the polymer already adsorbed. On the other hand, the "wiggling tails" of polymers 
can overlap. 

As in [3] there are two important quantities of interest: the number of covered mesh 
points A(t) and the number N(t) of polymers which are attached to the surface at time 
t. To understand the behaviour of A(t) and N(t), we introduce in Section[21an operator 
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Figure 2: One realization of pRSA algorithm for a — 0.8 and (3 — 0.1. The covered 
mesh points of the rectangular 100 x 100 mesh are shown at different times. 



formalism which makes it possible to derive a series expansion for N(t). We also derive 
series for A{t) and for numbers of point-polymers, cross-polymers and square-polymers 
adsorbed on the surface at time t. The operator formalism transforms the random 
sequential adsorption process into a set of combinatorial problems on the lattice. In 
some special cases, one can further simplify the resulting lattice combinatorial problems; 
we consider these special cases in Sectional The general problem is studied in Sectional 
To illustrate the precision of the derived formulas, we also provide a comparison of the 
results obtained by series expansion with those obtained by direct stochastic simulation, 
of particular interest is the time evolution of A(t) and N{t) and the dependence of the 
final adsorbed polymer layer on the parameters a and /?. We conclude with a discussion 
in Sectional 
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2 Operator formalism 



Let us denote by N{t) (resp. Np{t), Nc{t) and Ns{t)) the number of polymers (resp. 
point-polymers, cross-polymers and square-polymers) which are adsorbed on the surface 
at time t. Then we have 

Np{t) = (l-a- P)Nit), N,{t)=aN{t), Ns{t) ^ PN{t). (2.1) 

Let A{t) (resp. F{t)) be the number of covered (resp. vacant) mesh points at time t. 
Since dN/dt = F/AP and A = AP - F, we have 

Ait) = Al'(l-^it)]. (2.2) 



dt 

Let us define 

= lim N{t), N?° = lim iV„(t), = lim Nc{t), = lim 7V,(t). (2.3) 

Then p. 1(1 implies 

^ [l - a - I3)N°° , N^^aN°°, N^^PN°°. (2.4) 

Hence, the saturating values Nj^, and can be computed directly from N°°. 
Similarly, the time evolution of A{t), Np{t), Nc{t) and Ns{t) can be obtained from N{t) 
by H2.1|l - (|2.2|) . In this section, we develop an operator formalism framework to obtain 
the time evolution of N{t) and the limit N°°. Once we get N{t) and N°°, the rest of 
quantities of interest can be expressed by (|2.1|) , 1(2. 2|l and 1(2. 4|l and their dependence on 
the model parameters a and /3 can be also studied. 

In [20, an operator formalism was developed for studying the square lattice with 
nearest-neighbour exclusion. The results can be directly used to find an approximation 
of N°° for a — 1 and /3 = 0. If a = 1, then it is sufficient to keep track of the centers 
of cross-polymers. Each center of a cross-polymer excludes putting another center of 
a cross-polymer in the nearest neighbourhood of it. Hence, one can reformulate pRSA 
algorithm for a = 1 in terms of adsorption of points which excludes the nearest neigh- 
bourhood of them. Similarly, one can reformulate the pRSA algorithm as adsorption 
of points which excludes the nearest and the next nearest neighbourhood of them for 
a = and (3 — 1. However, if [a,/3] ^ {[0,1], [1,0], [0,0]}, then we have a mixture of 
polymers of different sizes in the solution and the approach of ^IS] cannot be directly 
used. In this section we present a generalization of the operator formalism for the case 
of arbitrary a and /3. 

We consider an M x M lattice (with periodic boundary conditions) to which polymers 
can adsorb. For each lattice point (i, j), we consider the state function G {0, 1, 2, 3}. 
Here, — means that lattice point is vacant or occupied by the "wiggling tail" 
of a cross-polymer/square-polymer (i.e. Vij = means that lattice point is free of 
centers of polymers/attached reactive groups), Vij = 1 means that lattice point («, j) is 
occupied by the point-polymer, Vij — 2 means that the lattice point is occupied by the 
center of the cross-polymer and i>ij = 3 means that the lattice point is occupied by the 
center of the square-polymer. Denoting 

|0) = [1,0,0,0]^ |1) = [0,1,0,0]^ |2> - [0,0,1,0]^ |3) = [0,0,0,1]^ 
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we identify every lattice point with the four- dimensional vector space M.^. Namely, the 
configuration of the M x M lattice will be expressed as 



M 

IK.})e0{|O),|l),|2),|3)}. 



The system state is given by 



(2.5) 



where the sum is taken over all possible configurations {I'ij} of the lattice and P{{i^ij}, t) 
is the probability of each configuration. It satisfies the normalization condition 



E P(K,},t) = i. 



(2.6) 



For each lattice point, we define cross, square and point annihilation operators 
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(2.7) 



More precisely, operator Ajj (resp. Bjj, and Cjj) acts as A (resp. B, and C) on 
the lattice point (i,j) and as the identity on all other lattice points. We also define 
creation operator A-^ (resp. bJ^-, and Ct^) as the transpose of operator Ajj (resp. 
Bi.j, and C^j). The cross-polymer number operator (resp. square-polymer number 
operator, and point-polymer number operator) is defined as = A\jAij (resp. 

N*j- = BtjBj^j, and Nf ^ = Cj^Cjj) which is, at the lattice point (^, j), a projection 
onto the one-dimensional subspaco spanned by vector |2) which corresponds to a cross. 
The "vacancy" number operator can be expressed as = A^ iAJ, = B^ iBJ — 



CijC\y Here, "vacancy" means that the lattice point is either free or covered by the 
tail of the cross-polymer /square-polymer, i.e. it is free of the attached reactive groups. 
We have N^.^- |*) 



)/6|^'). Let Ri j be an operator which 



is equal to the identity Id (resp. 0) operating on configurations in which lattice point 
(i, j) lies outside (within) the set of lattice points covered by tails of cross-polymers or 
square-polymers, i.e. 



°(Nr+i„-i 



■N? 



n: 
n; 



-i)(Nr- 



1 



-1 " 

N?-i,,-i), 



(2.8) 



where symbol o is used to emphasize that R^ j is the composition of operators which 
are typed on several lines (to simplify the resulting formulas, we skip the composition 
symbol o if the composed operators are typed on the same line). In the pRSA algorithm, 

we add the center of a cross-polymer (resp. square-polymer, and point-polymer) at the 
lattice site (i, j) at a rate a (resp. (3, and 1 — a — /?) if (i, j) is vacant and not covered by 
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the tail of a cross-polymer or square-polymer. This means that the state satisfies 
the master equation 



d_ 

dt 



(2.9) 



Solving ifT^ with the initial condition |{0}) = |0) ® |0) 

M 



|0), we obtain 



t 

M2 



^ (aAl_^ + /5Bt^ + (1 - a - /3)C^ - N1;^)r,„ 



|{0}). (2.10) 



Denoting \{u\) = |u) |u) ® • • • © |m) where \u) — [1, 1, 1, 1]-^, we can compute the 
number of polymers at time t by 



N(t) = m2({z,}|(N? , + , + N^,i)|*(t)). 



(2.11) 



Using (|2.10(l . we get 



7V(<) M^/ {u} 



(N?_i+N^^i + N^_i) 



o exp 



t 

M2 



M 

^ (aAl + /3Bt^. + (1 - a - /3)Cl^^. - N^,^.)r,, 



{0} 



{u} 



M 

J2 («Al_^. + + {l-a- (3)Cl - Nr,,)R.„ 



{0} 



1 / 1 



k=0 



kl \AP 



E n 



(aAt^.,,^+/?Bt^.,,^+(l-a-/?)Ct^,,^-N5;^.,,^.) 



R 



(2.12) 



{0}' 



Here, the last sum is done over all fc-tuples {{xj,yj)}j^i in the mesh. To evaluate this 
formula, let us note that we can consider the contributions of each mesh point separately. 
If (xi, yi) 7^ (1,1), then an operator of the following type acts on the mesh point {xi, yi): 



A^l.,y. + ^l.J'''m.,y. + ^l.,y. + N^.,,.]''^ ' ' ' ^ W,.,,., (2.13) 
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where 71, 72, 73, ■ • ■ , are nonnegative integers. Without loss of generality, we can assume 
73 > in what follows. The "building blocks" W^j-.j,. can be reasonably simplified if we 
take into account the following formulas: 

fN" +NP IfN" +NP 1 = N" +NP 

[N!^,,,, + NP,,,J[aAt_,,. + /3Bt^,,,+ (1 - a - P)Cl,yr ^l..y.\ = 



[^l.,y. + N^.y. + ^l.J [^^ky. + /^B^,,, + (1 - a - mkyr N^.yJ = 

If 71 = 72 = 0, then the building block (|2.13() can be rewritten in the form ±[q;AJ,. y. 
P^l^^y^+ (1 - a - P)Cl,,y,~'^l,,y^[^l,.xX' ■ We can easily observe that 



{u} 



{0} =0. 



Consequently, the first necessary condition for fc-tuple {(a^j, to have nonzero 

contribution to the formula H2.12|l is that for every {xi,yi) 7^ (1,1) in the fc-tuple, 
there must be j < z such that {xi,yi) is equal to {xj,yj) or one of its nearest or next 
nearest neighbours. In particular, we see that {xi, yi) = (1, 1) in order to have nonzero 
contribution of the fc-tuple { [xj , yj )}^=i- If 71 > 0, then (|TT^ satisfies 



{u} 



W 



Xi,yi 



{Q}) = {-ir{a + l3) 



(2.14) 



where we have denoted by 7^ the number of times that the mesh point (xi,yi) appears 
in the fc-tuple {{xj,yj)}j^i- Similarly, if 71 = and 72 > 0, then H2.13|l satisfies 



{u} 



W 



xi.yi 



{0} = (-1)^-/3. 



(2.15) 



Finally, considering the contribution of the first mesh point (xi, yi) — (1, 1), we get 



{■Nl, + Nl,+Nl,)[aAl,+PBl, + {l -a- /3)Cl,i -Nl'^J 



mj - (-i)''^-^ 

(2.16) 
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Let us define Vk as tlie set of all sequences s = {{xj,yj)}j^i, such that {xi,yi) = (1, 1) 
and for each i e {2, . . . , fc} there exists j < i such that {xi,yi) £ Sx^^yj, i.e. {xi,yi) 
is equal to (xj,yj) or one of its nearest or next nearest neighbours. Let us denote by 
uj{s) the number of distinct points in the sequence s € Vk- Let f(s) be the number of 
distinct points {xi,yi) e s, {xi,yi) ^ (1,1), satisfying that there exists j < i such that 
{xi,yi) € '^Xj.yj, i-e. ixi,yi) is equal to {xj,yj) or one of its nearest neighbours. Then 
we can rewrite H2.12|l (using (|2.14() - H2.16|l ) as 

k=i ' ^ ^ seVk 

Formula (|2.17|l is a starting point for the analysis of the pRSA algorithm. In order to 
evaluate coefficients of the series expansion (|2.17|l . we have to compute the quantities 

^ + for fc = 1,2,3,.... (2.18) 

Thus we have transformed the problem of the original pRSA algorithm to a combinatorial 
problem on the two-dimensional lattice. The problem can be further simplified if a = 
or /3 = as we will show in the following section. The general analysis of (|2.17|l for any 
a and /3 is given in Section 0] 



3 Analysis of pRSA algorithm in some special cases 



First, let us note that formula H2.17|l is consistent with the trivial case [a, /3] = [0, 0] 
We have 



fc=l 



t 

M2 



which is the exact formula ioi a — (3 ^ 0. This can be seen easily, since N(t) = 
A{t) = AP - F{t) where F{t) solves dF/dt = -F/AP. Formula l(TT7|l is also consistent 
with the cases = [1,0] and [a,/3] = [0,1] which were studied in QJIS]. Choosing 

[a,P] = [0,1], we get 



fc=i 



Vk 



(3.1) 



which is the formula derived in ^ Here, |7^fc| is the number of sequences in Vk- 
Similarly, if [a, (3] = [1,0], we get 



k=l 



Qk 



(3.2) 



where Qk is the set of all sequences s = {{xj,yj)}j^i, such that {xi,yi) = (1, 1) and for 
each i S {2, . . . , fc} there exists j < i such that {xi, yi) £ Xxj,yj, i-e. (xi, yi) is equal to 
{xj,yj) or one of its nearest neighbours. 
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The situation is more complicated if [a, f3] ^ {[0,0], [1,0], [0, 1]}. To evaluate coeffi- 
cients of the series H2.17|l . we have to compute the quantities ()2.18(l . If we use directly 
formula we would have to evaluate a different computationally intensive combi- 

natorial problem for each a and /3. Here we show that we can transform formula H2.17|l 
to the problem where computationally intensive part (involving Vk or Qk) is done in- 
dependently of a or f3. We start with the analysis of the pRSA algorithm in the special 
case a = 0. 



3.1 Special case a = 

If a = 0, then pRSA algorithm reduces to adsorption of point-polymers and square- 
polymers, and H2.17|l reads as follows 



k=l 



(3.3) 



If = 1, then (|3.3|l implies 1)3. It was observed in "5" that the Laplace transform can 
be used to further simplify the formula 1)3. Here, we show that the Laplace transform 
can help us to analyse H3.3I) for any f3. Taking the Laplace transform of (|3.3|) . term by 
term, we obtain 



N{u) 



iV(i)e~"*di 



Mix- 

11 ^—^ 



k=l 



(3.4) 



for sufhciently large u. Let us define Qk as the set of all sequences of k distinct points 
{(xj, yj)}*^^i, such that (xi,j/i) — (1, 1), and for each i g {2, . . . ,/c} there exists j < i 
such that {xi,yi) is equal to one of the nearest or the next nearest neighbours of {xj ,yj), 
i.e. {xi,yi) G Sooj,yj and {xi,yi) ^ {xj,yj). Then we have (using (|3.4() ) 



^ " = E ^r2 



k=l 



3 = 1 L 



iM2 



AP 

u 



k=l 



u ^ (fc-1)! 



dx. 



Taking the inverse Laplace transform, we obtain (for sufficiently large r) 



N{t) = 7^ N{r + m)e(''+™)*dM = AP 



l-cxp[-t/A/^] oo 



i-Px) 



k-1 



da;. (3.5) 



Let us define function 



^ix)^APJ2 



k=l 



i-pf-'x^' 

kl 



(3.6) 
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Then 1231) yields 



N{t) = * - exp 
In particular, the final coverage of the lattice can be computed as 



t 

"M2 



lim N{t) = lim l'(x) 



(3.7) 



(3. 



and other quantities of interest can be obtained by H2.ll) . (|2.2|1 and 1)2. 4|) . Thus, the 
adsorption algorithm has been reformulated to the problem of finding the numbers of 
sequences in the sets Qk, k = 1,2, . . . . Once, we have the numbers \Gk\ we can write ^'(x) 
for any /3 and compute N{t), Np{t), Ns{t) and A{t) by (ETjl, (E^J and provided 
that we can compute the sum of series (13 .61) with reasonable precision. To do so, we set 



(3.9) 



The first eight values of gk can be computed relatively easily as follows gi = 10000; 
£,2 = 40000; 53 = 146667; 54 = 500000; 55 = 1606667; ge = 4918889; 57 = 14461429 and 
58 = 41070290. Our task is to estimate the sum of series ()3.6|l knowing only the first 
eight partial sums 



5«(x) = E(-/3)'"' 



gkx 



= 1,2, 



(3.10) 



fc=i 



To do that, we use Shanks transformation computed by Wynn's algorithm |20[I19 | 
in the following way 



s„(x), for n — 1,2, 

l/{Sn{x) - Sn-l{x)) 



for n = 2, 3, . . . , 8, 



_n-fc+l 



-k-2 



{x) + l/{el-_r\^) 
for fc = 3, 4, 5, 6, 7, n = k, . 



£ 



n~k+l 
fe-1 



(^)) 



(3.11) 



To approximate the sum '^{x), we use the term £'j{x) which is also the Pade [4,3]- 
approximant since we use Shanks transformation for a power series |19) . Thus, we 
aproximate number of attached polymers as 



N{t)Ki £'^j{x{t)^, where x{t) 



1 — exp 



t 

"M2 



(3.12) 



The resuhs obtained by for /3 = 1 and M = 100 are given in Figure Ela) . To 

compute the time evolution of N{t) we chose an equidistant mesh for x in the interval 
(0, 1) and evaluated £7 by (|3.11() at each mesh point. Then the corresponding time t 
was computed by H3.12|l . In Figure|3Ja), we compare results obtained by approximation 
1)3.121) and by stochastic simulation of pRSA algorithm. We see that we get an excellent 
agreement between the theoretically derived formula and the simulation. The asymptotic 
coverage can be approximated as 



lim N{t) 

t — *oo 



e?(l) 



/34(1), and 7V„°°« (1-/3)4(1). 



(3.13) 
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Figure 3: pRSA algorithm for a — Q. (a) Time evolution of Nit) given by for 
a = 0, P = 1 and M — 100 (dashed line). Ten realizations of stochastic simulation of 
pRSA algorithm are plotted as thin solid lines, (b) iV°°, and as obtained by 

H5'. (solid lines) for AI = 100. We compare the approximate results with stochastic 
simulation of pRSA algorithm (20 realizations, each realization plotted as a dot). 



In Figure IHIb), we compare approximations of iV°°, N!^ and iV^ computed by H3.13|l 
with the results of stochastic simulations. We see that approximations (|3.13|) provide 
very good results for any j3. We can estimate relative error between approximation N'^^ 
and exact value iV°° as {N^^, - N°°)/N°°. Using we obtain N^^ - = 2.7 

and {N^p — N°°)/N°° = 0.15% for (3—1. Here, exact value of N°° was approximated 
by averaging over 100,000 realizations of the pRSA algorithm as N°° = 1869.8 for j3 ~ 1. 



3.2 Special case /3 = 

li f3 — 0, then the pRSA algorithm reduces to the adsorption of point-polymers and 
cross-polymers. The terms in the sum (|2.18|) are nonzero only if ^(s) — uj(s) — 1. Hence, 
H2.18|) can be rewritten as 

^ forfc = l,2,3,..., (3.14) 

where Qk is the set of all sequences s = {(a^j, %)}j=i, such that (xi, yi) = (1, 1) and for 
each i £ {2, . . . ,k} there exists j < i such that {xi,yi) G Xxj,yj- As before, uj{s) is the 
number of distinct points in the sequence s £ Qk. Using (|3.14|l , we can rewrite H2.17|l as 

oo / + \ *: 

m = E ^ (-1)'"' E (3-15) 

fe=i ' ^ ^ seQk 

Comparing formulas (|3.3ll and 13. 1511 . we find out only two differences: Vk in H3.3() is 
replaced by Qk in (|3.15() and /3 in l|3.3|l is replaced by a in H3.15|l . Consequently, taking 
the Laplace transform of 13.15(1 and using the same method as in Section ITTl we find 
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(compare with (|3.7|l ) 



N{t) = 1 - exp 



t 



where n{x) = ^ 

7, 1 



i-a) 



fc! 

k=l 



(3.16) 



where is the set of all sequences of k distinct points {(a^j, yj)}j=i, such that {xi,yi) = 
(1, 1), and for each i g {2, . . . , fc} there exists j < i such that {xi,yi) is equal to one 
of the nearest neighbours of (xj,yj), i.e. {xi,yi) € '^^j.j/j and {xi,yi) ^ {xj,yj). In 
particular, the final coverage of the lattice can be computed as 

N°° = lim N{t) = lim n{x) (3.17) 

and other quantities of interest can be obtained by H2.1|) . (|2.2|l and H2.4|) . Thus, the 
adsorption algorithm has been transformed to the problem of finding the numbers of 
sequences in the sets Tik, fc = 1, 2, . . . . Once, we have the numbers \Ti.k \ we can write D,{x) 
for any a and compute N{t), Np{t), Nc{t) and A{t) by l|TTB|) . ifTTIl and (|2!2l, provided 
that the series in fl{x) is convergent. It was pointed out in 6 that the convergence of 
series H3.1fi|l is slow for a = 1 and for x = 1. To overcome this difficulty, we could use 
Shanks transformation or Pade approximants as in Section f3.1l This approach works in 
general and we will use it in Section^ where the general analysis of pRSA algorithm is 
presented. Here, we present an alternative approach, rewriting series (I3.16|) in different 
variables. Several possibilities were shown and motivated in j^. Here, we write fl{x) as 

oo /"^ 3 
n{x{z)) =y^^akz'', where x= j -^—^^ ^r^d^. (3.18) 

Let us define 



1 1~L I 

hk^M^^-^. (3.19) 

To find hk, one has to solve a finite combinatorial problem. In this paper, we will 
make use of the first eight values of hk- They can be computed as follow hi — 10000; 
h2 = 20000; /i3 = 40000; = 73333; = 125333; = 202222; h^ = 311048; and 
h% = 459452. To find coefficients Ofc in H3.18|l . we substitute x = x{z) in H3.16|l . We 
differentiate the resulting series term by term eight times and we evaluate each derivative 
at z = to obtain: 

fli = hi 

02 = hi - a/i2 

03 = 2/11/3 — 2a/i2 + 

04 = hi/6-7ah2/3 + 3a^h3- a^hi (3.20) 
ag = -4/11/I5 - 5aft.2/3 + 5a^/i3 - 4a^/i4 + a^'/is 

ae = -4/11/9 -lla/i2/45 + lla^/i3/2-26a^/i4/3 + 5a''/i5-a^/i6 

07 ^ -20/ii/63 + 6a/i2/5 + 53a2/i3/15-38a3/i4/3 + 40a'*/i5/3-6a5/i6+a^^7 

as = 1.852 a/i2 - 13a^/i3/30- 63a^/i4/5 + 145a''/i5/6- 19a^/i6 + 7a^/i7 - a^/i8 
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Let z be a solution of equation 1 = 3 [1 + 2(1 — ^)''] ^d^ (one can numerically estimate 
z as 0.569). Moreover, let us denote 



(3.21) 



k=l 



where ai, . . . , og are given by H3.20|l . Then, using H3.18|l and (|3.16ll . we can approximate 
number of attached polymers as 



N{t) w Q[ z{t) ) , where z{t) is given by 1 — exp 



t 






-L 



1 + 2(1-0^ 



de 



(3.22) 

The results obtained by lliO^ for a = 1 and M = 100 are given in Figure ^a) . To 
compute time evolution of N{t), we chose an equidistant mesh in z-variable in interval 
[0,z] and evaluated fl by H3.21|l at each z. The corresponding time t was computed by 
H3.22I) . namely using the formula 



t = -M^ In 



1 - 



z{t) 



1 + 2(1-0^ 



In Figure 01a), we compare results obtained by approximation 13.22|l and by stochastic 
simulation of pRSA algorithm. We get a very good agreement between the theoretically 
derived formula and simulation. The asymptotic coverage can be approximated as 



iV°° = lim N{t) « n{z), w an{z), and w (1 - a) n{z). 



(3.23) 



In Figure Efb) , we compare approximations of N°° , and computed by H3.23|l 
with the results of stochastic simulations. We see that approximations (|3.23(l provide 
good results for any a. We can estimate relative error between approximation and 

= 5 and 



exact value N"^ as (N^^ - N°°)/N°°. Using ^^1^ . we obtain N^j, - N"" 
{N^p — N^)/N°° = 0.14% for a — 1. Here, N°° can be computed by averaging over 
many realizations of the pRSA algorithm as N°° = 3641 for a = 1. 

In this section, we used transformation of variables (|3.18|l to accelerate the conver- 
gence of series (I3.16f) . This transformation was suggested in ,6|| for pRSA algorithm 
with [a, /3] — [1, 0], but our analysis shows that it can give good results for any a. The 
problem with this approach in general is determining an appropriate change of variables. 
An easier, and more systematic, approach is to use a Shanks transformation or Fade 
approximants ^1 1201 ITT?) as we did in Section 13.11 and as we will do for the general 
analysis of the pRSA algorithm in Section 0] 



4 General analysis of pRSA algorithm 

To evaluate H2.17|l for general a and /?, we have to compute the quantities l|2.18|l for 
fc = 1,2,3,.... Direct evaluation of H2.18|l would require solving different combinatorial 
problems (weighted sums over all sequences in the set Vk) for different values of a and 
p. As in Section we show that a suitable reordering of terms can transform the 
set of combinatorial problems to only one combinatorial problem which can be solved 
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Figure 4: pRSA algorithm for 13 — 0. (a) Time evolution of N{t) given by ^3.2'^ for 
a — \, [3 — Q and M — 100 (dashed line). Ten realizations of stochastic simulation of 
pRSA algorithm are plotted as thin solid lines, (b) iV°°, and as obtained by 

^3. 23\i (solid lines) for AI = 100. We compare the approximate results with stochastic 
simulation of pRSA algorithm (20 realizations, each realization plotted as a dot). 



independently of the values of a and /3. Then the dependence of the number of attached 
polymers and number of covered binding sites can be easily studied. To do that, we first 
use the Laplace transform to rewrite (|2.17l) in terms of Gk- Here, as before Gk is the 
set of all sequences of k distinct points {{xj,yj)}j^i, such that {xi,yi) = (1, 1), and for 
each i e {2, . . . there exists j < i such that {xi,yi) g S^^^yj and {xi,yi) ^ {xj,yj). 
Following a similar analysis to that in Section ITTl we derive (compare with (|3.7|) ') 

(4.1) 
where 

CO / ^ 1 

<i>{x) ^ M^Y. ki + (4.2) 

k=l s^Qk 

where, as before, ^(s) is the number of distinct points {xi,yi) S s, {xi,yi) ^ (1,1), 
satisfying that there exists j < i such that {xi,yi) € Xxj,yj ■ Let (f>j, j — 1, . . . ,k, denote 
the number of sequences s G Gk, k — 1,2, ... , satisfying £,{s) = k — j. The numbers 0j 
for A: = 1, 2, . . . , 8, can be directly computed and they are given in Tabled Using the 
definition of (j>j, formula (|4.2|) can be rewritten to 

oc / -| \k — l ^ k 

Hx) = M' ^ ^~ \, E-^'t" + /3]'=-^/3^-^ (4.3) 
k=i ■ j=i 

Our task is to compute the sum of series (|4.3|) with reasonable precision, using only the 
first eight partial sums 

sr^ix) - J2 E '^'t" + f3f-'P'-\ n = 1, 2, . . . , 8. 

k=l ' j=l 



N{t) = $ 1 - exp 



t 
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j = 1 


j = 2 


j = 3 


j = 4 


j = 5 


i = 6 


j = 7 


i = 8 


k = 1 


1 
















k = 2 


4 


4 














k = 3 


24 


40 


24 












fe = 4 


176 


424 


424 


176 










k = 5 


1504 


4800 


6696 


4776 


1504 








fc = 6 


14560 


58368 


104752 


104280 


57640 


14560 






fc = 7 


156768 


761024 


1677680 


2135920 


1655336 


745064 


156768 




fc = 8 


1852512 


10603744 


27833952 


43206736 


42818768 


27137992 


10289192 


1852512 



(a) 



Table 1: Table of values of <j)j for k — 1,2, 

(b) 





Figure 5: (a) pRSA algorithm for M = 100 and a = 0.5; (b) pRSA algorithm for 
M = 100 and (3 = 0.5. In both cases, we present iV°°, , and as obtained by 
^4-51 (solid lines). We compare the approximate series expansion results with stochastic 
simulation of pRSA algorithm (20 realizations, each realization plotted as a dot). 



To do that, we use Shanks transformation computed by Wynn's algorithm H3.11|l and 
we approximate sum $(a;) by term £7(0;), as in Section im Thus we aproximate number 
of attached polymers as 



N{t) « e^ix{t) j , where x{t) is given by x{t) = 1 — exp 



t 

"M2 



(4.4) 



The asymptotic coverage can be approximated as 

N°°Ke^j{l), N^^ae^^il), 7V,°° « /3e?(l) and 7Vp°° w (1 a - /3) e?(l). (4.5) 

In Figure EJa), we compare approximations of N°° , N^, and computed by 
H4.5|l with the results of stochastic simulations for a = 0.5. The same plots for (3 = 0.5 
are given in Figure [S{b). We see that approximations (|4.5|l provide excellent results. 

The results obtained by H4.4|l for a = 0.5, /3 = 0.5 and M = 100 are given in Figure|H| 
To compute time evolution oi N(t), we chose an equidistant mesh in a;- variable in interval 
(0, 1) and evaluated £y by (|3.11|) at each x. Then the corresponding time t was computed 
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Figure 6: pRSA algorithm for a — 0.5, (3 — 0.5 and M — 100. (a) Time evolution 
of N{t) given by i4-4i (dashed line). Ten realizations of stochastic simulation of pRSA 
algorithm are plotted as thin solid lines, (b) Time evolution of A{t) given by ]4-4[ '^^'j^ 
Hi?.g|l (dashed line). Ten realizations of stochastic simulation of pRSA algorithm are 
plotted as thin solid lines. 



by (|4.4|l . To compute A{t) we used formula H2.2|l where the time derivative of N(t) was 
approximate by the backward-in-time finite difference of N(t). In Figure El we compare 
results obtained by approximation 14.4(1 and by stochastic simulation of pRSA algorithm. 
We see that we get a very good agreement between the theoretically derived formula 
and simulation. Finally, we present the time evolution of N{t) and A(t) for a = 0.8 
and P = 0.1 which is the situation shown in the illustrative computation in Figure |21 
In Figure [3 we compare results obtained by (|4.4|l with results obtained by stochastic 
simulation of pRSA algorithm. Again, we obtained an excellent agreement between the 
series expansion results and the stochastic simulation of the pRSA algorithm. 



5 Discussion 

In this paper we studied random sequential adsorption to the two-dimensional lattice. 
Our motivation was chemisorption from polydisperse solution of polymers. We gener- 
alized the operator formalism of j^E], derived series expansion results and presented 
efRcient methods to accelerate their convergence. In Section lTTl we used classical meth- 
ods for accelerating convergence of slowly converging series. In Section 11^.21 we also 
presented results obtained by a more specialized transformation of variables [S]. In both 
cases, the theoretical results compare well with the results of stochastic simulation of 
the pRSA algorithm. 

We assumed that the attached polymer can effectively shield a circle on the surface 
with radius r < 2h where h is the average distance between neighbouring binding sites. 
We worked with the rectangular mesh of binding sites to enable the reformulation of 
the problem in terms of the RSA on the rectangular lattice. One should view this 
simplification as a reasonable approximation of the problem where binding sites are 
more or less uniformly distributed on the surface. The restriction r < 2h can be also 
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Figure 7: pRSA algorithm for a — 0.8, (3 — 0.1 and M — 100. (a) Time evolution 
of N{t) given by i4-4i (dashed line). Ten realizations of stochastic simulation of pRSA 
algorithm are plotted as thin solid lines, (b) Time evolution of A{t) given by ]4-4[ '^^'j^ 
Hi?.g|l (dashed line). Ten realizations of stochastic simulation of pRSA algorithm are 
plotted as thin solid lines. 



relaxed and the operator formalism could be generalized to the case of a mixture of longer 
polymers too. However, one should have in mind that for larger r, the assumption that 
the "wiggling tails" of polymers can overlap has to be modified to take into account 
the higher probability to find the polymer chain close to the binding site; see for the 
general discussion of the polymer dynamics. 

Two-dimensional adsortption is more complicated to study because there is no sim- 
ple analogy of the exact approach which is available in one-dimension (see e.g. 0] or 
the integro-differential evolution equation framework which was used in [3]). More pre- 
cisely, one can formally write an evolution equation for the process (e.g. the master 
equation denoted (|2.9|l in this paper) but it can be solved only by various approxima- 
tion techniques For example, Nord et al ,12. study adsorption of dimers or larger 
connected sites of objects to two-dimensional lattice. They write a master equation in 
hierarchic form for conditional probabilities that a conditioned configuration of mesh 
points is empty given that some neighbouring conditioning sites are empty. Using a 
series of hierarchic truncation schemes |17| . they were able to estimate dynamics and 
saturating coverage of the adsorption process. The operator formalism presented is a 
useful alternative to methods based on approximate evolution equations. 

The theoretical treatment of irreversible polymer adsorption is given in |13| . They 
give a more detailed picture than is studied in this paper, by studying the structure of 
the resulting nonequilibrium layer in terms of the density profiles, and loop and contact 
fraction distributions. Adsorption of whole polymers to the surface, modelled as a self- 
avoiding random walk, was done in I18| where the results of Monte Carlo simulations 
are presented. It has been found that the coverage to its jamming limit is described 
by a power law t~^ where an exponent 7 depends on the chain length. In our case, 
we modelled the adsorption of polymers as adsorption of disks to the surface where 
the binding sites were arranged into the rectangular lattice. In particular, the presented 
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algorithm can be viewed as a generalization of the classical lattice RSA models. Random 
sequential adsorption has been subject of the intensive research for the last sixty years. 
The reader can find more details about the RSA in review articles 0] and 
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